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Abstract. In this paper we use new regularity and stability estimates for Alexandrov solutions to 
Monge- Ampere equations, recently established by De Philippis and Figalli [14] , to provide global in 
time existence of distributional solutions to the semigeostrophic equations on the 2-dimensional torus, 
under very mild assumptions on the initial data. A link with Lagrangian solutions is also discussed. 



The semigeostrophic equations are a simple model used in meteorology to describe large scale 
atmospheric flows. As explained for instance in [6j Section 2.2] and [18^ Section 1.1] (see also [11] for a 
more complete exposition), the semigeostrophic equations can be derived from the 3-d incompressible 
Euler equations, with Boussinesq and hydrostatic approximations, subject to a strong Coriolis force. 
Since for large scale atmospheric flows the Coriolis force dominates the advection term, the flow is 
mostly bi-dimensional. For this reason, the study of the semigeostrophic equations in 2-d or 3-d is 
pretty similar, and in order to simplify our presentation we focus here on the 2-dimentional periodic 
case, though we expect that our results could be extended to three dimensions. 

The semigeostrophic system on the 2-dimensional torus T 2 is given by 



and the functions ut and pt represent respectively the velocity and the pressure, while u\ is the so-called 
semi-geostrophic windQ Clearly the pressure is defined up to a (time-dependent) additive constant. In 
the sequel we are going to identify functions (and measures) defined on the torus T 2 with Z 2 -periodic 
functions defined on M 2 . 

Substituting the relation u\ = JVpt into the equation, the system (jl.ip can be rewritten as 



1. Introduction 




Here 



p is the initial datum, J is the rotation matrix given by 




(1.2) 



< 



' d t JVpt + JV 2 
V • ut = 

PQ = p° 



■p t u t + Vp t + Ju t = 



i 



Note that we are using the notation Ut, Pt, it? to denote the functions u(t, ■), p(t, ■), u 3 {t, ■) 
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with ut and p t periodic. 

Energetic considerations (see Section 3.2]) show that it is natural to assume that pt is (— 1)- 
convex, i.e., the function Pt(x) := pt(x) + |x| 2 /2 is convex on M 2 . If we denote with Jz? T 2 the (nor- 
malized) Lebesgue measure on the torus, then formalljQ pt := (V Pt)$Jfj2 satisfies the following dual 
problem (see the Appendix): 



( d tPt + v • (u tPt ) = o 

(13) \u t (x) = J(x-VP t *(x)) 

\Pt = (VP t )tt^ T2 

Here P t * is the convex conjugate of Pt, namely 

p t(y) : = su p (y- x - p t( x ))- 

Notice that, since Pt(x) — \x\ 2 /2 is periodic, 

(1.4) VP t (x + h) = VP t (x) + h VxgM 2 , h£Z 2 . 

Hence VP* can be viewed as a map from T 2 to T 2 and pt is a well defined measure on T 2 . One can 
also verify easily that the inverse map VP t * satisfies (|1.4|) as well. Accordingly, we shall understand 
(|l,3p as a PDE on T 2 , i.e., using test functions which are Z 2 -periodic in space. 

The dual problem f j 1 . 3 j) is nowadays pretty well understood. In particular, Benamou and Brenier 
proved in [6] existence of weak solutions to (|1.3p . see Theorem 13. li b elow . On the contrary, much less is 
known about the original system (jl.ip . Formally, given a solution pt of (|1.3p and defining P t * through 
the relation pt = (VP^jjJ^ 2 (namely the optimal transport map from pt to JZj2, see Theorem 12. lj) 
the pair (p t ,ut) given bjlij 



(1.5) 




P t (x) - \x\ 2 /2 

[dtVPt](VP t {x)) + [v 2 p;](vp 4 (x)) J(vp t (x) - x) 



solves (|1.2|) . However, being P 4 * just a convex function, a priori V 2 P 4 * is just a matrix- valued measure, 
thus as pointed out in |12j it is not clear the meaning to give to the previous equation. 

In this paper we prove that (11.5P is a well defined velocity field, and that the couple (pt, Ut) is a 
solution of (jl.ip in a distributional sense. In order to carry out our analysis, a fundamental tool is a 
recent result for solutions of the Monge-Ampere equation, proved by the third and fourth author in 
|14j . showing Llog fe L regularity on V 2 P 4 * (see Theorem 12. 2( ii) below). 

Thanks to this result, we can easily show that the second term appearing in the definition of the 
velocity Ut in fjl .5[) is a well defined L 1 function (see the proof of Theorem II. 2p . Moreover, following 
some ideas developed in [T7] we can show that the first term is also L , thus giving a meaning to 
ut (see Proposition I3.3|) . At this point we can prove that the pair (pt,ut) is actually a distributional 



Given a measure fi on T 2 and a Borel map / : T 2 -> T 2 , we define the measure /j/x through the relation 



Kv) d h»{y) = I h (f( x )) d ^( x ) 
x 2 Jt 2 



Because of the many compositions involved in this paper, we use the notation [dtf](g) (resp. [V/](p)) to denote the 
composition (dtf) ° g (resp. (V/) op), avoiding the ambiguous notation dtf(g) (resp. V/(<?)) 
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solution of system (|1.2p . Let us recall, following [12], the proper definition of weak Eulerian solution 
of (fl~2j) . 

Definition 1.1. Let p : T 2 x (0, oo) — > R and u : T 2 x (0, oo) — > R 2 . We say that (p, u) is a weafc 
Eulerian solution of (|1.2p if: 

- |u| G ^((C^oo),^^ 2 )), p G ^((O,^),^ 1 ' 00 ^ 2 )), and p t (x) + |x| 2 /2 is convex for any 
t > 0; 

- For every <j) G C C °°(T 2 x [0, oo)), it holds 

(1.6) J™ J JVp t (x)\d t <j) t {x) + u t (x) • V&(z)} 

- For a.e. i G (0, oo) it holds 

(1.7) / Vip(x) ■ u t {x) dx = 

We can now state our main result. 

Theorem 1.2. Let po : R 2 — > R a Z 2 -periodic function such that po(x) + |x| 2 /2 is convex, and 
assume that the measure (Id + Vpo)^ 2 is absolutely continuous with respect to S£ 2 with density po, 
namely 

(Id + Vpo)^ 2 = po^ 2 - 
Moreover, let us assume that both pq and 1/po belong to L°°(R 2 ). 

Let pt be the solution of (jl.3p given by Theorem I3.il and let Pt : R 2 — > R be the (unique up to 
an additive constant) convex function such that (VPt)^^ 2 = pt^£ 2 and Pt(x) — |x| 2 /2 is Z 2 '-periodic, 
P t * : R 2 — > R its convex conjugate. 

Then the couple (pt,ut) defined in (ll.5p is a weak Eulerian solution of (|1.2p . in the sense of Defi- 
nition \LT[ 

Although the vector field u provided by the previous theorem is only L , as explained in Section [5] 
we can associate to it a measure-preserving Lagrangian flow. In particular we recover (in the particular 
case of the 2-dimensional periodic setting) the result of Cullen and Feldman |12] on the existence of 
Lagrangian solutions to the semigeostrophic equations in physical space. 

The paper is structured as follows: in Section [2] we recall some preliminary results on optimal 
transport maps on the torus and their regularity. Then, in Section [3] we state the existence result of 
Benamou and Brenier for solutions to the dual problem (jl.3p , and we show some important regularity 
estimates on such solutions, which are used in Section 0] to prove Theorem 11.21 In Section [5] we prove 
the existence of a "Regular Lagrangian Flow" associated to the vector field u provided by Theorem 
11.21 Finally, in Section [6] we list some open problems. For completenes, in the Appendix we show the 
formal computation used to obtain (|1.3p from (II. 2p . 

Acknowledgement. L.A., G.D.P., and A.F. acknowledge the support of the ERC ADG GeMeThNES. 
A.F. was also supported by the NSF Grant DMS-0969962. 

2. Optimal transport maps on the torus and their regularity 

The following theorem can be found in [TO] (see for instance [141 Section 2] for the notion of 
Alexandrov solution of the Monge- Ampere equation) . 



^Vp t (x) + Jut{x]^4>t{x) dx dt 

+ I JVpo(x)(j)o(x) dx = 0; 
Jt 2 



for all if) G C°°(T ). 
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Theorem 2.1 (Existence of optimal maps on T 2 ). Let p and v be 1? -periodic Radon measures on M. 2 
such that p([0, l) 2 ) = v([0, l) 2 ) = 1 and p = pJ>? 2 with p > almost everywhere. Then there exists a 
unique (up to an additive constant) convex function P : M. 2 — > R such that (VP)^p = v and P — \x\ 2 /2 
is Z 2 -periodic. Moreover 

(2.1) VP{x + h) = VP(x) + h fora.e.x G R 2 , V/i G I?, 

(2.2) |VP(z) -x\< diam(T 2 ) = ^- for a.e. x G R 2 . 

In addition, if p = pS£ 2 , v = <rJ2? 2 , and there exist constants < A < A < oo such that A < p, a < A, 
then P is a strictly convex Alexandrov solution of 

det V 2 P{x) = f(x), with f(x^ - P ^ 



a{VP(x)) 



Proof. Existence of P follows from [10]. To prove uniqueness we observe that, under our assumption, 
also p*(y) := P*{y) — |y| 2 /2 is Z 2 -periodic. Hence, since 

P(x) = sup x • y — P*(y), 

y£R 2 

we get that the function p(x) := P(x) — |x| 2 /2 satisfies 

\y x\ 2 |yp 
p(x) = sup f - P*(y) + -|- 



\y + h — x\ 2 
sup sup ( p (y + h) 

ye[0,i\ 2 hez 2 



( d 2 T2 [x,y) \ 
sup ( p (y) I , 

yGT 2 v 



dj 2 (x,y) 



where dj2 is the quotient distance on the torus, and we used that p*(y) is Z 2 -periodic. This means 
that the function p is d 2 2 -convex, and that p* is its <i 2 2 -transform (see [20\ Chapter 5]). Hence 
VP = Id + Vp : T 2 — )■ T 2 is the unique (/x-a.e.) optimal transport map sending p onto v ( [191 
Theorem 9]), and since p > almost everywhere this uniquely characterizes P up to an additive 
constant. Finally, all the other properties of P follow from [10]. □ 

Combining the previous theorem and the known regularity results for strictly convex Alexandrov 
solutions of the Monge- Ampere equation (see [TJ [HJ [TO] HH [16] ) we have the following: 

Theorem 2.2 (Space regularity of optimal maps on T 2 ). Let p = p^ 2 , v = o-Jzf 2 , and let P be as in 
Theorem \2.1\ with J T2 P dx = 0. Then: 

(i) P € C 1,/3 (T 2 ) for some /3 = /3(A, A) G (0, 1), and there exists a constant C = C(A, A) such that 

\\P\\c^ < C 

(ii) P G W 2,1 (T 2 ), and for any k G N t/iere exists a constant C = C(X,A,k) such that 



[ \V 2 P\ log k + \V 2 P\ dx < C. 
Jt 2 
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(iii) If p,a 6 C fe - Q (T 2 ) for some k G N and a G (0,1), then P 6 C fe+2 ' Q (T 2 ) and i/iere exists a 
constant C = C(A, A, H^H^a, HcrH^.a) suc/i i/iat 

-P||cfc+2,a < C. 

Moreover, there exist two positive constants c\ and C2, depending only on A, A, ||p||c?o,a, and 
\\o\\(jo,a, such that 

Cl Id<X7 2 P(x) <c 2 Id VxGT 2 . 

3. The dual problem and the regularity of the velocity field 

In this section we recall some properties of solutions of (jl.3j) . and we show the L 1 integrability of 
the velocity field u t defined in (|1.5p . 

We know by Theorem 12.11 that pt uniquely defines Pt (and so also P£) through the relation 
(VPt)u-^r 2 = Pt up to an additive constant. We have the following result (see [6j [12]): 

Theorem 3.1 (Existence of solutions of (jl.3p ). Let Pq : IR 2 — > ]R be a convex function such that 
Po(x) — \x\ 2 /2 is 1? -periodic, (VPo)b-^t 2 ^ -^ 2 > an d the density po satisfies 0<A<po<A< 
oo. Then there exist convex functions P t , P t * : M 2 -)• R, with P t (x) - \x\ 2 /2 and P t *(y) - \y\ 2 /2 
periodic, uniquely determined up to time- dependent additive constants, such that (X7Pt)fiJ£ 2 = pt^£ 2 , 
(VP t *)j/9t = Ji?j2. In addition, setting Ut(x) = J(x — VP t *(x)), pt is a distributional solution to (jl.3p . 
namely 

(3.1) / / \ d t if t {x) + Vip t (x) ■ U t (x) \pt{x) dxdt + \ ip (x)p (x) dx = 

for every (p G C^°(1R 2 x [0,oo)) 1? -periodic in the space variable. 
Finally, the following regularity properties hold: 

(i) A < p t < A; 

(ii) p t J? 2 £C([0,oc),V w (T 2 ))B 

(iii) P t - f j2 P t , P t * - f j2 P t * G L°°([0, oo), J<' C °°(K 2 )) n C([0, oo), W^ r c (R 2 )) for every r G [1, oo); 

(iv) Moo < y/2/2. 

To be precise, in [61 [12] the proof is given in M 3 , but actually it can be rewritten verbatim on the 
2-dimensional torus, using the optimal transport maps provided by Theorem 12.11 Observe that, by 
Theorem 13. If ii). t h-> ptJz? 2 is weakly continuous, so pt is a well-defined function for every t > 0. 

Further regularity properties of VPt and VP 4 * with respect to time will be proved in Propositions l3.3l 
and 13.61 

In the proof of Theorem 11.21 we will need to test with functions which are merely W 1 ' 1 . This is 
made possible by the following lemma. 

Lemma 3.2. Let pt and Pt be as in Theorem\3Ji Then flU]) holds for every (p G W l ^(T 2 x [0,oo)) 
which is compactly supported in time. (Now <po(x) has to be understood in the sense of traces.) 

Proof. Let ip n G C°°(T 2 x [0,oo)) be strongly converging to </? in W 1 ' 1 , so that <Pq converges to tpo in 
L X (T 2 ). Taking into account that both pt and Ut are uniformly bounded from above in T 2 x [0, oo), 
we can apply (|3.ip to the test functions <p n and let n — > oo to obtain the same formula with (p. □ 



4 Here P W (T 2 ) is the space of probability measures on the torus endowed with the weak topology induced by the 
duality with C(T 2 ) 
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The following proposition, which provides the Sobolev regularity of 1 1— > VP/, is our main technical 
tool. Notice that, in order to prove Theorem II. 2\ only finiteness of the left hand side in (|3.2|) would be 
needed, and the proof of this fact involves only a smoothing argument, the regularity estimates of [H] 
collected in Theorem I2.2( ii). and the argument of \17\ Theorem 5.1]. However, the continuity result 
in [15] allows to show the validity of the natural a priori estimate on the left hand side in (I3.2p . 



Proposition 3.3 (Time regularity of optimal maps). Let pt and Pt be as in Theorem \3.1[ Then 
VP/ G W[ ' c (T 2 x [0,oo);]R 2 ), and for every k £ N there exists a constant C(k) such that, for almost 
every t > 0, 



(3.2) 



T 2 



p t \d t vp t *\\og k + (\d t vp t *\)dx 



< C(k) 



[ Pt \V 2 P t *\ log 2fe (|V 2 P;|)cfc + ess sup { Pt \U t \ 2 ) [ \V 2 P t *\dx). 

JT 2 T 2 JJ 2 / 



To prove Proposition 13.31 we need some preliminary results. 
Lemma 3.4. For every k EN we have 



- I 

e 



6 2 + 2 3(fc-i) a 2 log 2fc (a) V(a,6)G 



(3.3) a61og^(a6) < 2 fe - x 

Proof Prom the elementary inequalities 

Iog + (ta) < Iog + (t) + log + ( S ), (t + s) k < 2 k ~\t k + s k ) 
which hold for every t, s > 0, we infer 



log* (i) <[-) t 



ablog+(ab) < ab 



log. 



i k 



< 2 k ~ 1 ab 



< 2 



k-l 



< 2 



k-l 



+ 21og + (a) 



log^(-)+2 fc log^(a) 



^) 6 2 + 2 fc a61og^(a) 



6 2 + 6 2 + 22(fc -l) G 2 log 2 fc( 



□ 



which proves (I3.3p . 

Lemma 3.5 (Space-time regularity of transport). Let k 6 NU {0}, and let p S C°°(T 2 x [0, oo)) and 
U £ C°°(T 2 x [0,oo);M 2 ) satisfy 

< A < p t (x) < A < oo V (x, t) G T 2 x [0, oo), 

d tPt + V • (U tPt ) = mT 2 x[0,oo), 
and Jj 2 Ptdx = 1 for all t > 0. Let us consider convex conjugate maps Pt and P£ such that Pt(x) — 
\x\ 2 /2 and P t *(y) - \y\ 2 /2 are 1? -periodic, {VP?)§pt = (VP)^ 2 = Pt- Then: 
(i) P* t - f T2 P; G Lip loc ([0,oo);C fc (T 2 )) for any k G N. 
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(ii) The following linearized Monge-Ampere equation holds: 

(3.4) v • { Pt (v 2 p;r l d t vp;) = -v • ( Pt u t ). 

Proof. Let us fix T > 0. From the regularity theory for the Monge-Ampere equation (see Theorem l2.2j) 
we obtain that Pt G C°°(]R 2 ), uniformly for t G [0, T], and there exist universal constants ci, C2 > 
such that 

(3.5) aid < V 2 P t *{x) < c 2 Id V(s,t) G T 2 x [0,T]. 

Since VP 4 * is the inverse of VP t , by the smoothness of P t and (|3.5|) we deduce that P t * G C°°(]R 2 ), 
uniformly on [0, T]. 

Now, to prove (i), we need to investigate the time regularity of P 4 * — j-^Pt ■ Moreover, up to adding 
a time dependent constant to Pt, we can assume without loss of generality that jj 2 Pt = f° r au t- 
By the condition (VP*)$pt = Jzfjz we get that for any < s, t < T and x G R 2 it holds 

p s (x)- Pt (x) = det(V 2 P;(x))-det(V 2 P f *(x)) 
s — t s — t 



s - t 



Given a 2 x 2 matrix ^4 = (£jj)ij=i,2> we denote by M(A) the cofactor matrix of A. We recall that 

, ^ ddet(A) „ r , <N 

(3-7) / = Mij(A), 

and if A is invertible then M(A) satisfies the identity 

(3.8) M(A) = det(A) A~ l . 

Moreover, if A is symmetric and satisfies c\ld < A < c^ld for some positive constants ci, C2, then 

2 2 

(3.9) C -±Id < M(A) < ^Id. 

C2 Ci 



Hence, from ([321), dHTTJ) , (|33]) . and (pl9|) . it follows that 
(3.10) = £ (jf 1 A^(rV 2 P; + (1 - r)V 2 P t *) dr) % 



with 



c 2 f 1 c 2 

-^Jd < / Mm(tV 2 P.* + (1 - t)V 2 P;) dr < -2 Jd 

C2 Jo Ci 



Since V 2 P t * is smooth in space, uniformly on [0, T], by classical elliptic regularity theorjH it follows 
that for any k G N and a G (0, 1) there exists a constant C := C(||(p s — P t)/{s — i)||c fe > Q (T 2 x[o,T])) such 
that 

p;(x)-p t *(x) 



s - t 



< C. 

C k + 2 ' a (T 2 ) 



''Note that equation (|3.6|l is well denned on T 2 since p* — P s * is Z 2 -periodic. We also observe that P* — P s * has average 
zero on T 2 . 
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This proves point (i) in the statement. To prove the second part, we let s — > t in (|3.10p to obtain 

2 

(3.11) dtPt = Yl M ij(V 2 P?(x)) dtdijP t *(x). 

Taking into account the continuity equation and the well-known divergence-free property of the co- 
factor matrix 

J2d i M ij (y 2 P t *(x)) = 0, j = l,2, 

i 

we can rewrite (|3.11|) as 

2 

-V • (U tPt ) = £ diiMij^P^x^dtdjP^x)). 

Hence, using f|3.8j) and the Monge-Ampere equation det(V 2 P f *) = pt, we finally get (|3.4p . □ 

Proof of Proposition [373[ We closely follow the proof of \17\ Theorem 5.1], and we split the proof in 
two parts. In the first step we assume that 

(3.12) p t e C°°(T 2 x M), U t G C°°(T 2 x M; M 2 ) , 

(3.13) 0<A<p t <A<oo, 

(3.14) d t p t + V-(U t pt) = 0, 

(3.15) (VPf),^^^, 

and we prove that (|3.2p holds for every t > 0. In the second step we prove the general case through 
an approximation argument. 

Step 1: The regular case. Let us assume that the regularity assumptions (|3,12p . (|3.13p . (I3.14p . 
([3.15p hold. Moreover, up to adding a time dependent constant to Pt, we can assume without loss of 
generality that f j2 P t * = for alH > 0, so that by Lemma [33] we have d t P t * £ C°°(T 2 ). Fix t > 0. 
Multiplying (|3.4|) by dtP^ and integrating by parts, we get 



(3.16) 



/ pt\(v 2 p t *)- 1/2 dtVPt*\ 2 dx= [ Pt d t vp? ■{v 2 p;)- 1 d t vp? 

= - [ PtdtVPt* ■ U t dx. 
Jt 2 



dx 



(Since the matrix V 2 Pt*(x) is nonnegative, both its square root and the square root of its inverse are 
well-defined.) From Cauchy-Schwartz inequality it follows that the right-hand side of (|3.16p can be 
rewritten and estimated with 



(3.17) 



- f ptdtVPt* ■ (V 2 P t *r 1/2 (V 2 P t *)^ 2 U t dx 
Jt 2 

<([ p t \(V 2 P t *r 1/2 d t VP t *\ 2 dx] ([ p t \(V 2 P t *) 1/2 Ut\ 2 dx 
\Jt 2 / \Jt 2 

Moreover, the second factor in the right-hand side of (|3.17|) can be estimated with 

(3.18) I p t UfV 2 PtU t dx < max (pt\U t \ 2 ) [ \V 2 Pt\dx. 

Jt 2 t2 J T 2 



dx 
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Hence, from (pTTfi]) . (|3~TTjl . and (ETTKD it follows that 

(3.19) / p t \(V 2 P?y l l 2 d t VP?\ 2 dx<mB^{p t \U t \ 2 ) [ \V 2 P t *\dx. 

Jj2 T 2 Jf2 

We now apply Lemma E31 with a = \{V 2 P?) l / 2 \ and b = \(V 2 P t *)- 1/2 d t V 'P t *(x)\ to deduce the 
existence of a constant C(k) such that 

|a t vp;|iog^(|a t vp;i) < c(k) (\(v 2 p t *) 1/2 \ 2 iogf{\{v 2 p t *) 1/2 \ 2 ) + |(v 2 p t *)- 1/2 d t vp t *| 

= c(k) (|v 2 p;i io g ^(|v 2 p;i) + |(v 2 p;)- 1/2 5 t vp;i 2N 

Integrating the above inequality over T 2 and using ()3.19p . we finally obtain 

Pt \d t vp;\\o g k + {\d t vp;\)dx 

(3.20) <C(k)([ Pt \V 2 P t *\log 2 + k (\V 2 P t *\)dx+ [ p t \(V 2 P t T l/2 d t VP t *\ 

< C(k) ( [ Pt \V 2 P t *\ log^dV 2 ^!) dx + max { Pt \U t \ 2 ) [ |V 2 P t *| dx] , 
Vjt 2 tt 2 Jj2 ] 

which proves ()3.2[) . 

Step 2: The approximation argument. First of all, we extend the functions pt and Ut for t < by 
setting p t = po and Ut = for every t < 0. We notice that, with this definition, p t solves the continuity 
equation with velocity Ut on R 2 x R. 

Fix now (7i G C C °°(R 2 ), a 2 G C C °°(R), define the family of mollifiers (<r n ) neN as a n (x,t) := 
n 3 ai(nx)a2(nt), and set 

p * (7 

Since A < p < A then 

A < p n < A. 

Therefore both p n and C/ n are well defined and satisfy (|3.12p , (|3.13p , (I3.14|) . Moreover for every t > 
the function p™ is Z 2 -periodic and it is a probability density when restricted to (0, l) 2 (once again we 
are identifying periodic functions with functions defined on the torus). Let P™ be the only convex 
function such that (VP t n )jj,5? 2 = pf and its its convex conjugate P f n * satisfies J T2 P"* = for all 
t > 0. Since p" — > pt in L X (T 2 ) for any t > (recall that, by Theorem 13. If ii). pt is weakly continuous 
in time), from standard stability results for Alexandrov solutions of Monge- Ampere (see for instance 
[15]) it follows that 



(3.21) VP ( n * ^ VP ( * inL^T 2 ) 

for any t > 0. Moreover, by Theorems 12.11 and I2.2( ii). for every k G N there exists a constant 
C := C(A, A, fc) such that 

/ P ?\v 2 pr\iog k + (\v 2 pr\)dx<c, 

Jt 2 

and by the stability theorem in the Sobolev topology estabilished in [15\ Theorem 1.3] it follows that 

(3.22) / p?\V 2 Pr\log k + (\V 2 Pr\)dx^ [ pt\V 2 Pt*\logl(\V 2 P t *\)dx, 
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j2rm*\ j , / Iv72d* 



(3.23) / \V^P t n *\dx^ / \V^P t *\dx. 

Jt 2 Jt 2 

Finally, since the function (w,t) i— > F(w,t) = \w\ 2 /t is convex on M? x (0,oo), by Jensen inequality 
we get 

(3.24) l|p n |£H 2 ||oc = HF^tP.p^Hoo < ||p|C/| 2 ||oc. 

Let us fix T > and (p € C£°((0, T)) nonnegative. From the previous steps and Dunford-Pettis 
Theorem, it is clear that <j>{t) p" d t V 'P™* weakly converge to 4>(t) p t d t V Pf in iMY 2 x (0, T)). Moreover, 
since the function w i— > \w\\og+(\w\/r) is convex for every r G (0, oo) we can apply Ioffe lower 
semicontinuity theorem [H Theorem 5.8] to the functions <j)(t)pfdt'VPJ l * and <t>(t)pt to infer 

(3.25) 

/ ^(t) [ p t \d t VP;\\og k + {\d t VP;\)dxdt <liminf / <f>(t) f p^\d t VP^\\og k A\d t VP^\)dxdt. 
Jo Jt 2 n ^°° Jo Jj 2 

By Step 1 we can apply d£2J to p?,U?. Taking (I3T22D . (|5^5j> . ([3T21]) and (13351) into account, by 

Lebesgue dominated convergence theorem we obtain 

T (P(t) [ Pt \d t vp t *\iog k + (\d t vp?\)dxdt 



Jt 2 

<c(k) f T 4>(t) 

Jo 



[ /)i|V 2 P;|log 2 / c (|V 2 P;|)dx + esssup(p 4 |[/t| 2 ) / \V 2 P t *\dx] dt. 
J ■■• t 2 Jr 2 J 



Since this holds for every (p G C^°((0,T)) nonnegative, we obtain the desired result. □ 

It is clear from the proof of Proposition 13.31 that the particular coupling between the velocity field 
Ut and the transport map Pt is not used. Actually, using Theorem I2.2f ii) and |15} Theorem 1.3], 
and arguing again as in the proof of \17\ Theorem 5.1], the following more general statement holds 
(compare with [T71 Theorem 5.1, Equations (27) and (29)]): 

Proposition 3.6. Let p t and v% be such that 0<\<pt<A< oo, Vt 6 L^ C (T 2 x [0,oo),M 2 ), and 

d tP t + V • (v tPt ) = 0. 
Assume that Jj 2 Pt dx = 1 for all t > 0, let Pt be a convex function such that 

(VP t )^ T 2=p t J?j2, 

and denote by P t * its convex conjugate. 

Then VP t and VP t * belong to W£ c (T 2 x [0,oo);M 2 ). Moreover, for every k £ N there exists a 
constant C{k) such that, for almost every t > 0, 



(3.26) / p t \d t VP t *\\og k + (\d t VP;\)dx 
Jr 2 



< C(k) 



[ p t |V 2 P;|log^(|V 2 P;|)dx + esssup(p t |^| 2 ) / \V 2 P t *\dx 
Jt 2 t 2 Jt 2 

(3.27) / \d t VP t \log k + (\d t VP t \)dx 
Jt 2 

<C(k)( [ \V 2 P t \log 2 + k {\V 2 P t \)dx + esssu V (p t \v t \ 2 ) [ \V 2 P t 
\Jt 2 t 2 Jt 2 



2T3 t\dx 
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Proof. We just give a short sketch of the proof. Equation (I3.26f) can be proved following the same 
line of the proof of Proposition 13.31 To prove (|3.27p notice that by the approximation argument in 
the second step of the proof of Proposition 13.31 we can assume that the velocity and the density are 
smooth and hence, arguing as in Lemma 13. 5} we have that P u P t * £ Lip loc ([0, oo), C°°(T 2 )). Now, 
changing variables in the the left hand side of (I3.19P we get 

(3.28) I {[V 2 P;}(VP t )y 1/2 [d t VP t *](VP t ) 2 dx<mzx(p t \v t \ 2 ) f \V 2 P t *\dx. 

J T 2 T 2 Jj2 

Taking into account the identities 

[V^KVPi) = (V 2 ^)" 1 and [d t VP t *}(VP t ) + [V 2 P t *](V P t )d t V P t = 
which follow differentiating with respect to time and space VP t * o VP* = Id, Equation (|3.28p becomes 

[ \(V 2 P t y 1/2 d t VP t \ 2 dx<m&x( Pt \v t \ 2 ) [ \V 2 P t *\dx. 

J T 2 T 2 J T 2 

At this point the proof of (|3.27p is obtained arguing as in Proposition 13.31 □ 



4. Existence of an Eulerian solution 
In this section we prove Theorem 11.21 

Proof of Theorem 1 1.21 First of all notice that, thanks to Theorem I2.2f i) and Proposition 13.31 it holds 
|V 2 P;|, \d t VP t *\ G ^([Coo),^ 1 ^ 2 )). Moreover, since {VP t )^ T 2 = p t ^ T 2, it is immediate to check 
the function u in (fL~5l) is well-defmecfl and \u\ belongs to Lj^ c ([0, oo), L 1 (T 2 )). 

Let (f> S C^°(M? x [0, oo)) be a Z 2 -periodic function in space and let us consider the function 
ip : M 2 x [0, oo) — > IR 2 given by 

(4.i) My) ■= J(y - vP t *(y))MVP t *(y)). 

By Theorem 12.11 and the periodicity of <fi, pt(y) is Z 2 -periodic in the space variable. Moreover ipt 
is compactly supported in time, and Proposition 13.31 implies that p £ M /1,1 (]R 2 x [0, oo)). So, by 
Lemma 13.21 each component of the function <pt(y) is an admissible test function for (13. ip . For later 
use, we write down explicitly the derivatives of p: 



(4.2) 



d m (y) = -J[d t vp t *](y)M^Pt(y)) + J(v - vp t *( y ))[d t ^ t ](vp t *(y))+ 

+J( y - VPt(v))(\y<k](Pt(v)) • d t VP?(y)), 
y<pt(v) = J{Id - V 2 P?(y))MVP t *(y)) + J(y - Vi?(y)) ® ([V^](^*(y))V 2 J P,*(y)). 



Taking into account that (Vp^Jz^ = pt^f2 and that [V P^](V Pt(x)) = x almost everywhere, we 
can rewrite the boundary term in (|3.ip as 



(4.3) / (p (y)po(y)dy = / J(VP (x) - x)<f>o(x) dx = / JVp (x)<f>o(x) dx. 

Jt 2 Jt 2 Jb? 



^Note that the composition of V 2 P* with VPt makes sense. Indeed, by the conditions (VPt)tt^r 2 = ptJ^i 2 "C -2r 2 i 
if we change the value of V 2 P t * in a set of measure zero, also [V 2 P t *](VPt) will change only on a set of measure zero. 



12 L. AMBROSIO, M. COLOMBO, G. DE PHILIPPIS, AND A. FIGALLI 

In the same way, since Ut{y) = J(y — VP/(y)), we can use (14. 2D to rewrite the other term as 
(4 ' 4 i 

J J J2 {dm(y) + v<pt(y) ■ u t (y)} Pt (y) d y dt 

= I" I {- J[dtVP t *](VP t {x))(l> t (x) + J(VP t (x) - x)8tMx) 
Jo Jt 2 l 

+ J(VP t (x) - x)(V&(a;) • [dtVP?](VP t (x))) 

+ [J(ld - V 2 P;{VP t {x)))<f> t {x) + J(VP t (x) - x) ® (V T 4> t {x)V 2 P*(VP t {x)))] J(VP t (x) - x)} dxdt 



which, taking into account the formula (|1.5p for u, after rearranging the terms turns out to be equal 
to 



(4.5) / / {JVp t {x){d t <t) t {x) + u t {x)-Vct>t{x)) + {-Vp t {x)- Ju t {x))(t>t{x)} dxdt. 
Jo Jt 2 

Hence, combining (|4.3p . (|4.4p . (|4.5p . and (|3.ip . we obtain the validity of fjl .6j> . 

Now we prove (jl.7p . Given G C^°(0,oo) and a Z 2 -periodic function ^ G C°°(R 2 ), let us consider 
the function ip : R 2 x [0, oo) -> M defined by 

(4.6) ^(y) := 0(t)^(VP t *(y)). 

As in the previous case, we have that 93 is Z -periodic in the space variable and ip G W 1 ' 1 ^ 2 x [0,oo)), 
so we can use tp as a test function in ()1.7p . Then, identities analogous to (j4.2j) yield 

/"OO p 

0= / {dtiptM + VtptM-UtWyptWdydt 
Jo Jt 2 

r 

4>'{t) / ip(x)dxdt 
Jt 2 

+ j <f>(t) / |v^(x) • 9 t VP t *(VP t (x)) + V T V(x)V 2 P;(VP(x))J(VP t (x) - x)} dxdt 

<p(t) / X7ip(x) ■ u t (x) dx dt. 
Jt 2 

Since (ft is arbitrary we obtain 

/ Vip(x) ■ ut(x) dx = for a.e. t > 0. 

By a standard density argument it follows that the above equation holds outside a negligible set of 
times independent of the test function ip, thus proving (|1.7p . □ 

5. Existence of a Regular Lagrangian Flow for the semigeostrophic velocity field 

We start with the definition of Regular Lagrangian Flow for a given vector field b, inspired by [21 [3]: 

Definition 5.1. Given a Borel, locally integrable vector field b : T 2 x (0, 00) — > M 2 , we say that a 
Borel function F : T 2 x [0, 00) — > T 2 is a Regular Lagrangian Flow (in short RLF) associated to b if 
the following two conditions are satisfied. 

(a) For almost every x G T 2 the map 1 1-> -Fj(x) is locally absolutely continuous in [0, 00) and 

(5.1) F t {x)=x+ I b s {F s {x))dx Wg[0,oo). 

Jo 
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(b) For every t G [0, oo) it holds (P t )#Jzf T 2 < CJrf T 2, with C G [0, oo) independent of t. 
A particular class of RLFs is the collection of the measure-preserving ones, where (b) is strengthened 

to 

(F t ) # J? T 2 = if T 2 Vt > 0. 

Notice that a priori the above definition depends on the choice of the representative of b in the 
Lebesgue equivalence class, since modifications of b in Lebesgue negligible sets could destroy condition 

(a) . However, a simple argument based on Fubini's theorem shows that the combination of (a) and 

(b) is indeed invariant (see [21 Section 6]): in other words, if b = b a.e. in T 2 x (0,oo), then every RLF 
associated to b is also a RLF associated to b. 

We show existence of a measure-preserving RLF associated to the vector field u defined by 

(5.2) u t (x) = [d t VP t *](VP t (x)) + [V 2 P 4 *](VP(x))J(VP(x) - x), 

where Pt and P t * are as in Theorem ll.2l Recall also that, under these assumptions, \u\ G L^ c ([0, oo), L 1 (T 2 

Existence for weaker notion of Lagrangian flow of the semigeostrophic equations was proved by 
Cullen and Feldman, see \12\ Definition 2.4], but since at that time the results of |14j were not 
available the velocity could not be defined, not even as a function. Hence, they had to adopt a more 
indirect definition. We shall prove indeed that their flow is a flow according to Definition 15.11 We 
discuss the uniqueness issue in the last section. 

Theorem 5.2. Let us assume that the hypotheses of Theorem \1.2\ are satisfied, and let Pt and P£ be 

the convex functions such that 

Then, for Ut given by (|5.2p there exists a measure-preserving RLF F associated to Uf. Moreover F 
is invertible in the sense that for all t > there exist Borel maps F£ such that F^(F t ) = Id and 
F t (F t *) = Id a.e. in T 2 . 

Proof. Let us consider the velocity field in the dual variables Ut(x) = J(x — VP t *(x)). Since P t * 
is convex, U t G BV(T 2 ;R 2 ) uniformly in time (actually, by Theorem E^ii) U t G H /1,1 (T 2 ; ]R 2 )). 
Moreover Ut is divergence-free. Hence, by the theory of Regular Lagrangian Flows associated to BV 
vector fields PIS], there exists a uniqu^ measure-preserving RLF G : T 2 x [0, oo) -> T 2 associated to 
U. 

We now define^ 

(5.3) F t (y):=VP t *(G t (VP (y))). 

The validity of property (b) in Definition 15 . 1 1 and the invertibility of F follow from the same arguments 
of [121 Propositions 2.14 and 2.17]. Hence we only have to show that property (a) in Definition 15.11 
holds. 

Let us define Q n := B * a 11 , where B is a Sobolev and uniformly continuous extension of VP* 
to T 2 x R, and o~ n is a standard family of mollifiers in T 2 x R. It is well known that Q n — > VP* 



^ The uniqueness of Regular Lagrangian Flows has to be understood in the following way: if Gi, G2 : T 2 x [0, 00) — > T 2 
are two RLFs associated to U, then the integral curves Gi(-, x) and G2(-, x) are equal for «5f 2 -a.e. x. 

^Observe that the definition of F makes sense. Indeed, by Theorem I2.2f i). both maps VPo and VP t * are Holder 
continuous in space. Morever, by the weak continuity in time of t pt (Theorem I3.1f ii^ and the stability results for 
Alexandrov solutions of Monge- Ampere, VP* is continuous both in space and time. Finally, since (V-Fb)n-2r2 <K J&i?, if 
we change the value of G in a set of measure zero, also F will change only on a set of measure zero. 
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locally uniformly and in the strong topology of (T 2 x [0, oo)). Thus, using the measure-preserving 
property of Gt, for all T > we get 

0= lim f ^ \\Q2 -VP*\ + \d t Q? -d t VP*\ + \VQ™ -V 2 P*\\dydt. 

n^ooJ T2 J Q L J 

lim / / {|Q?(G t )-Vi?(G t )| + |[$Q?](Gy^^ 



n— >oo J^2 

Up to a (not re-labeled) subsequence the previous convergence is pointwise in space, namely, for almost 
every x G T 2 , 

[ T {\QnG t (x)) - VP t *(G t (x))\ + \[d t Q?](G t (x)) - [d t V P t *](G t (x))\ 
(5.4) Jo 1 

+ |[VQ?](G t (x)) - [V 2 P t *](G t (x))\}dt -> 0. 

Hence, since G is a RLF and by assumption 

(VP )~S? T 2 «^ T 2, 

for almost every y we have that (|5.4p holds at x = VPo(?/)> an d the function i i— > Gt(x) is absolutely 
continuous on [0, T], with derivative given by 

jG t {x) = U t {G t {x)) = J(G t (x) - VP t *(G t (x))) for a.e. t G [0,T]. 

Let us fix such an y. Since Q n is smooth, the function Q™(Gt(x)) is absolutely continuous in [0,T] 
and its time derivative is given by 

±(Q?(G t (x))) = [d t Q?](G t (x)) + \yQ?](G t (x))J(G t (x) - VP t *(G t (*))). 

Hence, since J(Gt(x) — VP t *(Gi(x))) = U(Gt{x)) is uniformly bounded, from (j5.4j) we get 
(5.5) 

lim ±(Q?( Gt (x))) = [d t VP t *](G t (x)) + [V 2 P t *](G t (x))J(G t (x) - VP t *(G t (x))) := v t (y) in L l (0,T). 
Recalling that 

lim QUG t (x)) = VP?(G t (x)) = F t (y) Vt € [0,T], 

n— >oo 

we infer that Ft{y) is absolutely continuous in [0,T] (being the limit in H^ 1 ' 1 (0,T) of absolutely 
continuous maps). Moreover, by taking the limit as n — > oo in the identity 

Qt(G t (x)) = QUG (x)) + J* ±{Q^{G T {x))) dr, 

thanks to (15.51) we get 

(5-6) F t (y) =F (y)+ [ v T (y) dr. 

Jo 

To obtain (|5.ip we only need to show that vt(y) = ut{Ft{y)), which follows at once from (|5.2p . (|5.3p . 
and ([531) . 

□ 
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6. Open problems 

In this short section we point out some open problems. The first one is of course uniqueness for the 
Cauchy problem, both at the level of (jl.3j) and at the level of (|l,2p . Let us point out that a priori the 
two problems are not equivalent, because we proved that solutions to (jl.3p induce solutions to (II. 2p . 
but at the moment the converse implication is only formal (see the Appendix). 

Another open question is the uniqueness of the regular Lagrangian flow associated to u. Uniqueness 
is known, thanks to the results in [2], for the flow G in the dual variables with velocity Ut(y) = 
J(y — VP t *(y)); actually, in light of the Llog k L Sobolev regularity of U, even the quantitative stability 
results of [13] are by now available for G. We were able in the previous section to prove that flows 
Gt of U induce flows Ft of u, via the transformation Ft = VP ( * oG t o VPo. However, our proof used 
the boundedness of U, an information we do not have when we try to reverse the implication, namely 
that regular Lagrangian flows F of u induce regular Lagrangian flows G of U via the transformation 
G = VPt o Ft o VPq . This question could be settled, at least in the class of measure- preserving 
Lagrangian flows, if the following conjecture had a positive answer: 

Conjecture. Let f G W^ l {(0,T) x T 2 ;M 2 ) n C([0,T] x T 2 ;]R 2 ), and let H t be a measure-preserving 
Lagrangian Bow relative to b. Assume that 

(6.1) \dtft]{H t (x)) + [V ft)(H t (x))bt{Ht(x)) G L\0,T) for a.e. x G T 2 . 

Then for a.e. x G T 2 the map t h-> ft{H t {x)) is absolutely continuous. 

In our case, / = VP and Ht is a measure-preserving flow associated to b = u; with these choices, 
the term in (|6.ip is equal to Ut(x), so it is even bounded, even though the summands in the expression 
might be unbounded. 

We remark that if we assume that / G W 1,q for some q > 1, and that 



f f ^-H t {x) dtdx= f I \b t {x)\ p dtdx < oo, 
Jt 2 Jo dt J T 2 J 



P 



1' 



then a simple approximation argument based on convolving / with smooth convolution kernels, as the 
one used in the proof of Theorem 15.21 provides a positive answer to the above conjecture. (This result 
can also be seen as a particular case of the general theory of weak gradients and absolute continuity 
along curves recently developed in [HE]. However, if / is not continuous, one needs to replace / with 
a suitable "precise representative" in its Lebesgue equivalence class.) Observe that, in this latter case, 
(|6.ip is automatically satisfied by Young inequality. 

Appendix A. From physical to dual variables 

For completeness, we formally show how the dual equation (jl.3p is derived from system (jl.2p . Taking 
into account the definition of Pt, the identities J 2 = —Id, Vpt(y)+y = VPt(y), V 2 pt(y)+Id = V 2 Pt(y) 
and the fact that ut is divergence- free, for every test function <p we obtain 

d f d f f d 

- <p(x)dpt(x) = - v{VP t (y))dy= V<p(VP t (y))- -Vp t (y)dy 
dt J T 2 at Jj2 Jj2 dt 

V<p(VP t (y)) ■ {(V 2 Pt (y) + Id)u t (y) - JV Pt (y)} dy 



f V[p(VPt{y))]-u t (y)dy+ [ V^(VP t (y)) • J(VP t (y) - y) dy 

I Vif(x) ■ J(x- VP t *(x))d Pt (x) = I V(p(x)-U t (x)dpt(x). 
Jt 2 Jt 2 
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Notice that this formal derivation holds independently of u (only the divergence-free condition of u 
is needed), and that u does not appear explicitly in (jl.3p . 
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